Recursive partitioned inversion of large (1500 x 1500) symmetric matrices

A recursive algorithm was designed to invert large, dense, symmetric, positive definite matrices using small amounts of computer core, i.e., a small fraction of the core needed to store the complete matrix. The described algorithm is a generalized Gaussian elimination technique. Other algorithms are also discussed for the Cholesky decomposition and step inversion techniques. The purpose of the inversion algorithm is to solve large linear systems of normal equations generated by working geodetic problems. The algorithm was incorporated into a computer program called SOLVE. In the past the SOLVE program has been used in obtaining solutions published as the Goddard earth models

Improvement in the geopotential derived from satellite and surface data (GEM 7 and 8)

A refinement was obtained in the earth's gravitational field using satellite and surface data. In addition to a more complete treatment of data previously employed on 27 satellites, the new satellite solution (Goddard Earth Model 7) includes 64,000 laser measurements taken on 7 satellites during the international satellite geodesy experiment (ISAGEX) program. The GEM 7, containing 400 harmonic terms, is complete through degree and order 16. The companion solution GEM 8 combines the same satellite data as in GEM 7 with surface gravimetry over 39% of the earth. The GEM 8 is complete to degree and order 25. Extensive tests on data independent of the solution show that the undulation of the geoidal surface computed by GEM 7 has an accuracy of about 3m (rms). The overall accuracy of the geoid estimated by GEM 8 is estimated to be about 4-1/4m (rms), an improvement of almost 1m over previous solutions

Goddard earth models (5 and 6)

A comprehensive earth model has been developed that consists of two complementary gravitational fields and center-of-mass locations for 134 tracking stations on the earth's surface. One gravitational field is derived solely from satellite tracking data. This data on 27 satellite orbits is the most extensive used for such a solution. A second solution uses this data with 13,400 simultaneous events from satellite camera observations and surface gravimetric anomalies. The satellite-only solution as a whole is accurate to about 4.5 milligals as judged by the surface gravity data. The majority of the station coordinates are accurate to better than 10 meters as judged by independent results from geodetic surveys and by Doppler tracking of both distant space probes and near earth orbits

Gravitational field models for the earth (GEM 1 and 2)

Two models of the earth's gravitational field have been computed at Goddard Space Flight Center. The first, Goddard Earth Model 1 (GEM 1), has been derived from satellite tracking data. The second, Goddard Earth Model 2 (GEM 2), has been derived from a combination of satellite tracking and surface gravimetric data. The geopotential models are represented in spherical harmonics complete to degree and order 16 for the combined solution and complete to degree and order 12 for the satellite solution. Both solutions include zonal terms to degree 21 and related satellite resonant coefficients to degree 22. The satellite data consisted primarily of optical data processed on 300 weekly orbital arcs for 25 close earth satellites. Surface gravity data were employed in the form of 5 deg x 5 deg mean free-air gravity anomalies providing about 70% world coverage. Station locations were obtained for 46 tracking sites by combining electronic, laser, and additional optical tracking data with the above satellite data. Analysis of the radial positions of these stations and a value of mean gravity on the geoid indicated a mean equatorial radius for the earth of about 6378145 meters. Results of geopotential tests on satellite data not used in the solution show that better agreement was obtained with the GEM 1 and GEM 2 models than with the 1969 Smithsonian Standard Earth 2 model